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Abstract 

The scattering of a Dirac particle has been studied for a quaternionic potential step. In the 
potential region an additional diffusion solution is obtained. The quaternionic solution which 
generalizes the complex one presents an amplification of the reflection and transmission rates. A 
detailed analysis of the quaternionic spinorial velocities sheed new light on the additional solution. 

For pure quaternionic potentials, the interesting and surprising result of total transmission is 
found. This suggests that the presence of pure quaternionic potentials cannot be seen by analyzing 
the reflection or transmission rates. It has be observed by measuring the mean value of some 
operator. 

I. INTRODUCTION 

There are several proposals to formulate quantum mechanics m by generalizing the complex number 
field. The most obvious generalization is done by introducing quaternions [3H2], but there are other 
possibilities, like Clifford algebras [S] and p—adic numbers [5]. Even inside each proposal, there may 
be internal divisions. In this article, we study a quaternionic quantum mechanics endowed with a 
quaternionic scalar product [5]. 

Our interest in the quaternionic formulation of Dirac equation is motivated by the recent interesting 
results obtained in the standard Dirac theory for the diffusion [T0l[Tl], tunneling [12] and Klein 
zones [13]. The generalization to quaternionic potentials |7] could be useful to understand many 
hidden aspects of the Dirac solutions. Even simple solutions, as those involving the potential step, 
contain interesting results. In the case of quaternionic quantum mechanics, there are fewer studies 
which investigate the quaternionic formulation of Dirac equation. In Adler’s book [5], the theoretical 
foundations have been laid down, nevertheless explicit physical examples are still quite rare. More 
recently, the Dirac solutions in the presence of quaternionic potentials have been detailed discussed 
in ref. |7]. In the potential region, the quaternionic solutions allow two values for the wave-function 
momentum. This represents the first important difference with respect to the complex case in which 
only one value is found. A better understanding of the additional quaternionic solution can be 
achieved by studying an explicit example of quaternionic potential and calculating the reflection and 
transmission coefficients. In order to do this, we recall the general solution presented in a previous 
article [7] and study the scattering by a quaternionic potential step. In the complex scattering, for 
one-dimensional motion, we have not spin flip in the reflected and transmitted beams pifTI]. Thus, 
the spin flip is not considered too for the quaternionic case. The additional quaternionic solutions and 
the new reflection and transmission coefficients satisfy conservation of probability and in the complex 
limit reproduce the standard Dirac scattering. A novel and intriguing behavior is observed for the 
scattering by a purely quaternionic step. In this case, there is no reflection and the wave-function 
shows an oscillating spatial pattern similar to the one observed for circularly polarized light. 

The article is organized as follows. In section II, for the convenience of the reader we fix our 
notation and recall the Dirac solution in the presence of a quaternionic potential. In section III, we 
obtain the reflection and transmission coefficients and calculate the probability current. In section 
IV, we discuss the effect of the quaternionic potential on the particle velocities. The last sections 
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round off the article with a brief discussion on the total transmission for pure quaternionic potentials 
(section V), final considerations and future possible investigations (section VI). 


II. THE QUATERNIONIC DIRAC EQUATION 

As usual in quaternionic quantum mechanics [3], the Dirac equation is written in terms of an anti- 
hermitian Hamiltonian operator. In this spirit, we introduce the anti-hermitian quaternionic potential 
step 

{0 for z < 0 (region I) , z Vo + j + fc ^2 for z > 0 (region II) } 
where are real constants. In the potential region (region II), the Dirac equation then reads 

9t'P(r,t) = -\a-V + imP + iVo+ +kV 2 ] (I) 


For Vi = 14 = 0, the complex step is recovered. The Dirac Hamiltonian conatins the matrices a and 
P which satisfy the algebra 


a. = a.' 


/3 = ^t 


oP = = 1, {/3,q;} = 0 


and 


{ar, Or} = 2 Srs 1- (2) 


In the Dirac representation [2], we have 


a = 


0 cr 

(T 0 


/3 = 


1 0 
0 -1 


(3) 


where cr are the well known complex Pauli matrices. For convenience of notation, the quaternionic 
potential can be rewritten as follows 

Wq = V 2 + iVi = \Wo\ exp [i arctan(I4/t4) ] = |IUo| exp [z </>] . 

Considering a motion along the z-axis, p = (0 , 0 , Q ), introducing the the quaternionic spinor 

= [u + jw] exp[i{Q z - Et)] , (4) 


where u and v are complex spinors, we obtain two coupled complex matrix equations, 

{E — Q — m P — Vo) u = —WgW 

{E — Qa^ + mP +Vo) w = —WoU . (5) 


The first equation is obtained from the complex part of equation o, and the second one is obtained 
taking the pure quaternionic part. Equations ([5]) can be rewritten as eigenvalue equations, 

{E-Qa3 + mP+Vo){E-Qao-mP-Vo)u = |M4| u, 

[E -Qa^-mp - Vo)(E -Qao + mP +Vo) w = |H4| w. (6) 


The determinants of the matrices on the left hand side of equations ([6]) are equal, and consequently 
the eigenvalues obtained from any of them are also equal. Non-trivial solutions have to satisfy 

r 21 


det 


{E - Q as + m P + I 4 ) {E -Qa^-mP - Vo) - |M4 


= 0 . 


(7) 


From equation 0, we obtain two squared momenta for particles moving in the region H, namely 


Ql=ql+ |IUor± 2(5, 


( 8 ) 


where 

= {E ±Vof - m^, 5 = ^ {EVof + (p|M4| f - EVo and p = \/E"^ - wP . 
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From equation (jS]), > 0, and Q may be either positive or negative. Hence, we find three energy 

zones determined by the Q_ value, 


1) diffusion zone 

q[ >0 

and 

E > 

VI Wo 

^ + {Vg+my , 


2) tunneling zone 

1 

A 

o 

and 

max 

m, yj 

\Wg\^ +[Vg-mf 

<E<yJ\Wg\^ + {Vg + mY 

3) Klein zone 

> 0 

and 

E < 

^\Wg 

to 

+ 

A 

1 

to 


These regions generalize the diffusion, tunneling, 

[IQHI3], 

and Klein zones obtained for the complex case 

1) 

diffusion zone 

2 

<7_ 

> 0 

and 

E > Vg + m , 


2) 

tunneling zone 

2 

< 0 

and 

Vg — m<E<Vo + m, 

3) 

Klein zone 

2 

> 0 

and 

E <Vo — m . 



For diffusion nnmi] and Klein m zones, we find oscillatory solutions interpreted as particles or 
anti-particles. In the tunneling zone |12] the wave function becomes evanescent and thus we have not 
particle or anti-particle propagation. 

In order to explicitly determine the quaternionic wave functions, we calculate the eigenvectors of 
equations The four component eigenvectors can be conveniently written by defining the two- 
dimensional vectors x = {(!> 0)*) (0) !)*}■ After simple algebraic manipulation, for Q = Q_, we 
find 

-1 

and w = — Wo {E — Q_ + m[3 + K,) u , (9) 


and u = —Wo (E — ao — m j3 — Vq) w . (10) 


u = 




and for Q = , 

Q+ 


E + Vo+ m+ 


0-3 X 


0-3 X 


III. THE MATCHING CONDITIONS 


In order to study the scattering of a Dirac particle by the quaternionic step, for simplicity of calculation 
we choose a spin orientation for the incoming particle, i.e. X = (1 0)* explicit in the complex 
spinors u, w, u and w their dependence on Q, Vg and Wg- For example, in region I, we have an incident 
particle with momentum p = — m^, thus represented hy u[Q = p; Vg = 0 , Wg = 0]. Obviously, 

in the same region there are reflected particles with momenta —p, thus represented by m[— p; 0,0]. 
There is no spin flip in the reflection or transmission, thus we can eliminate the zero component in 
our wave functions and work with two-component quaternionic spinors. 

In region I (potential free region), the quaternionic wave function is then given by 


'I'j(z) = u[p; 0, 0] exp[ip 2 :]-I- u[—p;0,0]R + jw[—p;0,0]R exp[—ipz], 


( 11 ) 


where R and R are complex coefficients. Note that the presence of complex exponential in (jT]) requires 
multiplication from the right by complex coefficients. In region II (potential region), the quaternionic 
wave function is characterized by two possible momenta, <5+ s-nd Q-^ 3'iid its explicit expression is 


= u[Q_; Vg, Wg] +j w[Q_; Vg, Wg] T exp[iQ_ z] + 

u[<3+; Vo, Wg] +jw[Q+-, Vg, Wg] T exp[j(5^ z] . 


( 12 ) 
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As happens for the reflection coefficients, the transmission coefficients T and T are also complex. The 
matching condition, 'I'j(O) = d>^j(0), implies 

u[p- 0, 0]+u[-p; 0, 0] i? = u[Q_; K, ITo] T + ^Q^; K, ITo] f 

w[-p-, 0,0] R = w[Q_-, Fo, Vho] T + w[Q^-, TTo] f . (13) 

From these coupled equations, we can obtain the complex coefficients R, R, T, T in terms of the wave 
function parameters and consequently calculate the transmission and reflection rates. 


III.A THE REFLECTION AND TRANSMISSION COEFFICIENTS 

To simplify our notation let us introduce the adimensional quantity 

a = 


P 


E + m 

and explicitly rewrite the quaternionic wave function in region I, 


= 

and region II, 


exp[ipz] + 


1 

— a 


R exp[ — ipz] + j 


where 

= 


1 

A 


-jWo 


M_ 

N 


T eyip[iQ_z] + 




Qi 


E ± Vo + m ± 


, M.= 


Q±A^ + E - TO =F K 


M. 




R exp[ — ipz] , 


A. 


T exp[iQ^z] , 
(E + TO Vq)A^ + 


ql-Ql 


The continuity equations at z = 0 lead to the following system 


1 

—a 

—a 

1 


R = -Wo 


T -Wo 


M. 


M_ 

N 


T 


T 


T . 


(14) 


R = 

T = 


After algebraic manipulations, we find 

(fl — A_ )(a + A^ ) + ]Wo I (AT^ — a ) {M_ + a N_ ) 

(a + A_)(a + A^)-|ITor(M^ +aA^)(M_ +aN_) ’ 

2 Q, (a A^ ) 

(a + AJ(a + A^)-|Wor(M^ +aA^)(M_ +aN_) ’ 

M -A^N 

^ ^ ^ ^ 

CL + 

M +aN 

T = Wo ~ T . 

a + 

In the complex limit. Wo —)> 0, we recover the standard reflection and transmission rates [2], 

Rc = Tc = 0 , 


(15) 


(16) 


and 


Rc 

Tr 


p {E - Vo + m) - {E + to) 
p (E - W + to) + (E + to) 
2p(E -Vo + m) 
p (E - Vo + to) + g__ (E + to) 


(17) 
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III.B THE CURRENT PROBABILITY DENSITY 


The relation between the reflection and transmission coefficients can be obtained by using the conti¬ 
nuity equation 

| + V.J = 0, (18) 

where 

p — and J — . 

Observing that ti'iti' is independent of time and the spatial dependence of dt is only on the z-coordinate, 
we find 

^'■''03 ^' = 0 . 

The continuity of the current density between the free region I and the potential region II implies 

+ K a, 4 -^ + a, 4 -^ = a, 03 • 


For the previous equation, we find 

^ A_ + |VFo|'m*Y_ 


|R| +\R\ 


2 a 


+ h.c. \T\ + 


2 f A^ + lITol 


2 a 


-bh.c. |T| = 1 . (19) 


In the complex limit, IFo —>■ 0, we recover the well know continuity equation 

(U + q*JiE + m) 


\Rc 


2p {E -Vo + m) 


\Tc\ =1. 


( 20 ) 


For E > Vo-|-m, we have gM. and pc > 0. Eq. (EOl) then implies a reflection rate < 1 (diffusion zone). 
For Vo — m < E < Vq + m, have q_ < 0. Consequently Pc = 0 and Ea. (|20)) implies a reflection rate 
= 1. Finally, for E < Vo — m, we have q_ € R and pc < 0, we find a reflection rate > 1 suggesting pair 
production mM- In Fig. 1, we plot the reflection rates for different value of quaternionic potentials 
as a function of E/m. In the Klein zone, the reflection rate increases by increasing the quaternionic 
potential. In this energy zone, the quaternionic perturbation thus contributes to the phenomenon of 
pair production. This is confirmed by the transmission rates in Fig. 2. It is important to observe that 
for complex tunneling only evanescent waves, exp[—|g_|z], appears in the potential region and this 
is seen by observing (Fig.la) that the reflection rate is 1. The presence of a quaternionic potential 
breaks total reflection and this is explained by the fact that together evanescent wave, exp[—1(5_ |z], 
also oscillatory waves, exp[i|(5+l-z], are present in the potential region. The numerical calculation for 
the combined reflection and transmission coefficients is shown in Fig. 3 and confirms the continuity 
equation IH]). From Fig. 3, we also observe a shift in the starting point of the diffusion zone and more 
important we have a decreasing energy zone for the evanescent waves by increasing the quaternionic 
potential. The last observation clearly suggest that oscillatory waves kill evanescent waves in presence 
of great quaternionic potentials. Finally, we also find that the reflection and transmission rates |i?| 
and |r| are very small compared to the reflection rates |i?| and |r|, thus as a first approximation we 
can use these last coefficient rates to study quaternionic perturbations of complex quantum scattering 
problems. 


IV. GROUP VELOCITIES ANALYSIS 


The incoming complex plane wave has the standard oscillatory exponential, exp[i{pz — Et], and 
consequently its group velocity is given by 


dE p 
dp E 


( 21 ) 
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The transmitted quaternionic plane waves contain, in the diffusion and Klein zones, the oscillatory 
exponentials 

exp [i {Q^z - Et)] , 

which determine the following group velocities 


Recalling that 


dE 

d(5± 


qI = {E±vS -m^ + \wS ±2 + IIToH - m-|fTor - EV, ^ 


we immediately obtain 


Q±‘^Q± 


E±Vo± 


E{v",+\Wo\") 



dE , 


and consequently 



1± 


K +|W"ol 




( 22 ) 


The additional quaternionic spinor is then characterized by a velocity which is a typical velocity 
of diffusion, see Fig. 4a. Indeed, we have not tunneling or Klein zones as happens for the velocity v _, 
see Fig. 4b. By increasing the quaternionic part of the potential we see decreasing the velocity of the 
additional quaternionic spinor and by increasing the incoming energy we obtain an increasing value 
of the velocity. These are typical effects of diffusion. As observed before, the velocity v_ presents 
three different energy zones (see Fig. 4b), the diffusion zone for E > \/|IF|o + {Vq + m) where the 
potential acts decelerating the particles, the Klein zone for E < \/|IF|o + (fo ~ ’m) where the negative 
velocities solutions generates a flux from the right to the left which gives a reflection rates greater than 
one and consequently interpreted in terms of pair production and finally the intermediate tunneling 
zone where we have not propagation due to the presence of evanescent solutions. It is interesting to 
observe that in the diffusion zone, we have two propagation regimes and consequently two different 
velocities. These different velocities permit to decouple the two quaternionic solutions. In the pure 
quaternionic limit (V), —>-0), as we shall discuss in our conclusions, these velocities tend to the same 
value and this value is p/E the free propagation velocity. 


V. TOTAL TRANSMISSION 


Our analysis shows a surprising effect for pure quaternionic potentials, Vq = 0. In this case. 


Q^=p ± |IFo|, = a, 


M.=t 


|VFo 


and 


N^=T 


1 


and the continuity equation leads to the following system for the reflection and transmission coefh- 
cients, _ _ _ 

l±Rq=Tg+Tg and ±Rg = -Tq + Tg . (23) 

The solutions are immediately obtained by showing total transmission, 

Rg=Rg=0 Tg = ^ and fg=TgE^. (24) 

The transmitted wave-function for the pure quaternionic potential is then given by 


4- 


q,T 


cosdWok) — k 


sin(|IFo|z) 


^ipz 


cos(|Wok) sin(IIFok) 

-kE* sin(\Wo\z) cosdlFoi^) 


(25) 
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and propagates with group velocity p/E. Due to the fact that the transmitted quaternionic spinor 
satisfies 

^ = 1 , 
q,T q,T ’ 

and moves with the same velocity of the free incoming particle seems impossible to see the presence 
of pure quaternionic potentials. Nevertheless the problem of an apparent invisible presence of pure 
quaternionic potential can be solved by taking the mean value of operators. For example, the z 
component spin mean value for the complex wave function is given by 


^ . 2 l + a" 2(l + a") 2F; ’ 

Cjinc Cjinc V / 


(26) 


which in the non-relativistic limit, p << to, reproduces the well-known result (S^)^ = 1/2. In the 
presence of a pure quaternionic potential, we should find 




S 'S . 

g,inc * q,inc 


^ (1 a) 


cos(2|Wo|z) — sin(2 |Wo| z) 
sin(2|Woiz) -cos(2iWoiz) 


2 1 + a^ 

|Wo| z) = ^ cos(2 |Wo| z) , 


(27) 


which implies a z-dependence on the mean value of spin similar to what happens for circular polarized 
light [14]. 


VI. CONCLUSIONS 

In this article we have studied the scattering of a Dirac particle by a quaternionic potential step. Our 
analysis shows the presence of an additional quaternionic solution characterized by a group diffusion 
velocity . The quaternionic solution which generalizes the standard one is characterized by a group 
velocity v_ and presents three energy zones, diffusion-tunneling-Klein zones. For diffusion the two 
solution are, in general, different and could be decouple the two solution in the potential region 
showing the presence of a quaternionic potential. By increasing the quaternionic part of the potantial 
we increase the Klein zone and decrease the tunneling zone, see Fig. 4b. 

The most interesting effect is seen for a pure quaternionic potential where total transmission and 
a z dependence on the spin operators appear. This suggests a possible investigation in quaternionic 
optics [6]. Before to do it, due to the fact that we have to introduce Gaussian laser, in a forthcoming, 
we aim to extend the plane wave analysis done in this article to quaternionic wave packets |5|. 
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Figure 1: The reflection rates as a function of the incoming energy for a fixed complex potential Vg = 
3m (blue continuous line) and for quaternionic perturbation \Wq\ = (l,2,3)m (cyan/magenta/red 
dashed and dotted lines). By increasing the quaternionic part of the potential, the tunneling zone 
decreases. We also observe that in presence of quaternionic perturbation, due to the additional 
diffusion solution, we have not total reflcetion in the tunneling zone. 
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Figure 2: The transmission rates are plotted as function of the incoming energy for a fixed complex 
potential and different quaternionic perturbation. As observed in the caption of Fig. 1, in presence of 
quaternionic potentials total reflection breaks down. This is clear from the plot in (b) which confirms 
the presence of additional diffusion solutions. 
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REFLECTION AND TRANSMISSION 



Figure 3: The reflection and transmission rates satisfy the probability conservation. In the Klein zone, 
we find a negative transmission flux and consequently a reflection rate greater than one. This effect 
is explaind in terms of pair production. The quaternionic perturbation amplify the pair production. 
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DIRAC PARTICLE VELOCITY 
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Fig. 4 


Figure 4: The group velocities, plotted as function of the incoming energy, confirms the analysis done 
for the reflection and transmission rates. The quaternionic perturbation acts in the diffusion zone by 
decelerating the particle and by decreasing the tunneling energy zone. 
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